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discharge  (Vg  - 1.3  KeV,  Ig  - 63  ma)  at  low  pressure  (P  - 2-4(10"6)Torr)  and  low  magnetic 
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orders  of  magnitude  larger  than  would  be  expected  for  a crow-field  coUisional  diffusion  process, 
a result  which  is  identified  with  the  turbulent  state  of  collective  beam-plasma  interactions. 


PLASMA  DIFFUSION  IN  A 


SP ACE- S I MULAT ION  BE AM-PLASM A-D I SCHARGE 

I.  INTRODUCTION 

Beam-plasma  interactions  have  received  increased 
attention  in  recent  years  largely  as  a result  of  consid- 
erations bearing  on  vehicle  neutralization  during  space- 
borne  accelerator  experiments,  enhanced  beam-plasma 
ionization  processes,  and  in  general,  collective  phenomena 
initiated  by  beam  injection  into  neutral  gas  and  charged- 
particle  environments.  With  plans  for  Shuttle-borne 
experiments  directed  at  controlled  beam-plasma  interactions 
and  the  use  of  beams  to  probe  ionospheric/magnetospheric 
electric  and  magnetic  fields,  it  became  important  to 
conduct  laboratory  simulations  which  supported  the  Shuttle 
plans  and  helped  establish  definitive  Shuttle-borne  experi- 
ments. 

As  part  of  just  such  a plan  involving  a continuing 
investigation  of  large-facility  beam-plasma  interactions 
(Bernstein,  et  al.  1975,1977,1978)  the  Naval  Research 
Laboratory  (NRL)  participated  in  a recent  series  of  experi- 
ments at  the  NASA  Johnson  Space  Center  (JSC).  The  NRL 


Note:  Manuscript  submitted  Dtctmbtr  11,  1978. 


contribution  involved  the  direct  measurement  of  electron 
density  and  temperature  under  varying  conditions  of  beam- 
plasma  and  neutral  gas  parameters.  Reported  here  are 
the  first  large  facility  profiles  of  electron  density  and 
temperature  under  conditions  of  enhanced  beam-plasma 
ionization ...  conditions  which  have  come  to  be  known  as  the  BPD, 
the  "beam-plasma-discharge"  (Getty  and  Smullin,  1963; 

Bernstein  et  al.,1978).  The  observed  density  profile  is 
found  to  be  in  good  agreement  with  a two-dimensional  diffu- 
sion model  and  this  result  makes  possible  the  first  estimate 
for  the  cross-field  electron  diffusion  coefficient  in  a 
large  facility  BPD.  In  the  following  Sections  the  experi- 
ment is  described,  the  theoretical  model  is  defined,  and 
the  diffusion  coefficient  is  derived. 

II.  EXPERIMENTAL  RESULTS 

The  Experiment.  The  technique  of  the  pulsed-plasma- 
3 

probe  (P  ) was  employed  to  measure  the  beam-plasma  profile. 

3 

The  P provides  simultaneous  measurements  of  N , T , V , 
and  <5N  (-*"P  (k) ) and  is  particularly  useful  under  dynamic 
plasma  conditions  and  in  environments  that  can  contaminate 
electrode  surfaces  (Holmes  and  Szuszczewicz,  1975; 
Szuszczewicz  and  Holmes,  1975,1976,1977).  Both  these 
conditions  prevailed  to  various  degrees  in  the  JSC  experi- 
ments. 

The  experimental  configuration  (similar  to  that 
described  by  Bernstein  et  al  1978)  is  illustrated  in 
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Figure  1.  The  P was  mounted  on  a traversal  mechanism 
positioned  at  approximately  8 m above  the  injection  point  of 
the  beam.  The  tungsten  cathode  gun  was  operated  at  63  ma 
and  1.3  KV  DC  with  the  chamber  pressure  in  the  2-4(10-6) 
torr  range.  A combination  of  coil  current  and  the  Earth's 
field  established  vertical  and  horizontal  component  B- 
fields  at  Bz=1.21  gauss  and  Bx_y=>0.27  gauss,  respectively. 
Under  these  conditions  a BPD  was  established  and  the  probe 
traversal  mechanism  was  exercised  to  determine  the  plasma's 
radial  profile  out  to  4.6  m.  The  resulting  measurements  of 
T and  N are  presented  in  Figure  2. 

Following  the  diffusion-model  density-profile 
predictions  (discussed  below)  the  experimental  curves 
were  fit  with  exponential  distributions  found  to  be 

Te  [lO3  °k]=  6.93  exp  (-r/1.99)  (1) 

Ne  ClO6  cm-3  ] « 7.01  exp  (-r/2.49)  (2) 


III.  THEORETICAL  CONSIDERATIONS 
For  purposes  of  comparison,  a simple  two-dimensional 
diffusion  model  is  developed  which  follows  the  approach  used 
in  the  study  of  small  laboratory  arc  plasmas  by  Bohm  et  al.  , 
(1949a).  The  approach  suggests  the  consideration  of  a 
differential  plasma  element  of  thickness  dr  and  of  length  h 
parallel  and  equal  to  the  length  of  the  beam  as  depicted 
in  Figure  3.  The  differential  element  is  taken  outside  the 
domain  of  the  primary  beam  and  is  assumed  to  be  homogeneous 
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along  its  length.  The  equilibrium  density  then  arises  from 
the  balance  of  charge  accretion  by  differential  drift  across 


the  element,  i.e.,  j (r)  - j (r+dr),  and  end  losses  (Z=0 
and  h)  due  to  axial  diffusion.  Requiring  dNe/dt  = 0 and 
assuming  an  axial  diffusion  loss  directly  proportional  to 
density  (i.e.,  differential  end  losses  * 2 aeN  dr),  the 

balance  between  production  and  loss  requires 

(dje/dr)  h<>r»  -2a  l Ne6  r (3) 

which  is  rewritten  as 

dje/dr=- 2 a®  Ne/h  . (4a) 

Similar  arguments  for  the  plasma's  ion  component  yield 

dj./dr  = -2  a*  N±/h  (4b) 


Independent  of  Eqs.  (3)  and  (4)  the  diffusion  currents  Je 
and  j ^ can  be  described  by  their  respective  mobility 


equations 
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where  V is  an  electric  potential  and  D and  D.  are  the 

6 X 

perpendicular  electron  and  ion  diffusion  coefficients, 
respectively.  Taking  the  one-dimensional  space-derivative 
d/dr  of  Eqs.  ( 5)  and  substituting  (4)  yields 
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when  Ne  = and  the  terms  have  been  rearranged.  Addition 

of  Eqs.  (6)  eliminates  the  potential-dependent  term  and 
results  in 


= 2 h 
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with  a solution 
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e I e 
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N°  exp  (/SR)  - N* 
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In  the  limit  R <*>  , the  solution  becomes 


where 


Ne  - N”  exp  (-0r) 

^ , 2(a zDiTe  - tXTi) 
DeDi(Ti  + Te>  h 


(10) 


IV.  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT 
The  agreement  between  theoretical  prediction  (Eq.(9)) 
and  experimental  results  is  found  to  be  good  (see  bottom 
panel,  Fig.  2)  and  provides  the  opportunity  for  determining 
the  effective  electron  diffusion  perpendicular  to  the  magnetic 
field,  i.e.,  the  radial  diffusion  in  the  model  of  Fig.  3. 

To  accomplish  this,  consider 


0-1  »/Deh/o  eV  * 2.49 
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where  the  value  on  the  right  side  is  taken  from  the  exponential 

in  Eq.  (2).  In  addition,  consider  as  the  effective 
velocity  of  axially  diffusing  electrons  which  have  sufficient 
energy  to  overcome  chamber  sheath  potentials  and  be  com- 
pletely lost  to  the  system.  In  this  case, 
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where  is  the  chamber-sheath  potential.  With  Vc  ■ 

( 4 . 8+? * ? ) kT_/e  (Chen,  1965;  Szuszczewicz.  1972)  and 

•1  • w 0 "" 

T - 3 (103)  °K  Eq.  (12)  yields 
e 

a ®(nom.)  - 7.0  (104)  cm  /sec  (13a) 

z 

a ®(max)  ■ 2.8  (105)  cm  /sec  (13b) 

z 
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oi  (min)  ■ 4.3  (10  ) cm  /sec  (13c) 
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Substituting 

the 

z 

values  into  Eq.  (11)  provides 

D^(nom) 

V 

N 

to 

• 

to 

(106) 

cm2 /sec 

(14a) 

D (max) 
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00 

00 
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cm  /sec 

(14b) 

D^(min) 

- 1.3 

<106) 

cm2 /sec 

(14c) 

as  the  first  experimental  estimate  of  the  effective  electron 
diffusion  coefficient  perpendicular  to  the  magnetic  field 
in  a BPD. 

The  values  for  are  orders  of  magnitude  larger  than 
would  be  expected  for  cross-field  collisional  diffusion. 

This  result  in  itself  is  not  surprising  since  a variety 
of  observations  indicate  that  the  diffusion  of  plasma 
across  magnetic  lines  of  force  can  be  substantially  enhanced 
above  the  collisional  diffusion  rate  when  fluctuating 


fields  are  present  in  the  plasma  (Spitzer . 1962;  Kadomtsev, 
1965).  This  is  indeed  the  case  in  a BPD  which  has  as  one  of 


its  characteristics  the  presence  of  fluctuating  electric 
fields  that  are  associated  with  large  amplitude  plasma 
waves.  The  enhancement  of  diffusion  in  turbulent  plasmas 
goes  back  to  Bohm  et  al.  (1949b)  who  semiempirically  de- 
scribed the  rate  by  what  has  come  to  be  known  as  the  Bohm 
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where  £t  j = e V and  [B]  = gauss.  For  our  nominal 

temperature  of  3(103)  °K  and  a total  field  B = 1.2  gauss 

z 

2 

Eq.  (13)  yields  Dg  = 1.3(10  )cm  /sec.  This  calculation 
has  not  been  intended  to  confirm  or  deny  the  validity 
of  Eqs.  (14)  but  has  been  presented  as  the  only  possible 
comparison  with  an  existing  concept  of  turbulent  diffusion. 
While  the  Bohm  formula  appears  to  apply  to  a surprising 
number  of  different  experiments,  care  must  be  exercised 
so  that  the  illusion  of  universal  validity  does  not  auto- 
matically develop.  It  is  evident  that  a first-principles 
derivation  of  the  coefficient  of  turbulent  diffusion  can- 
not be  obtained  without  a detailed  investigation  of  the 
plasma,  particularly  the  fluctuation  power  spectra  and 
associated  instability  processes  (Kadomtsev.  1965; 
Papadopoulos , private  communication  1978).  The  experi- 
mental and  theoretical  aspects  of  these  details,  as  they 
relate  to  the  JSC  experiments,  are  currently  underway 
at  NRL  (Szuszczewicz  and  Papadopoulos,  private  communi- 
cation 1978). 
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SUMMARY 


The  first  experimental  radial  profiles  of  electron 
density  and  temperature  in  a large  facility  beam-plasma- 
discharge  are  found  to  fit  exponential  functions  defined  by 


Te  [°K] 


6.93  (1CT)  exp  (-r[  m]  /1.99) 
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NJcm-3]  = 7.01  (106)  exp  (-r[  m ] /2.49). 

While  the  present  stage  of  research  does  not  provide  a 
theoretical  basis  for  the  temperature  distribution,  the 
electron  density  results  do  agree  with  a two  dimensional 
diffusion  model  which  predicts 

Ne  ■ Ne°  exp  | -(2«®/D^h)*r  } 

in  the  limiting  case  of  large  chambers  with  T./T  — *0.  With 

I C 

taken  as  the  effective  velocity  of  axially  diffusing 
z 

4 

electrons  (assigned  a nominal  value  of  7 (10  ) cm/sec)  the 
identity  of  the  experimental  profile  with  the  theoretical 
prediction  yields 

DJ’e(nom)  = 2.2(103)  cm3 /sec 

for  the  radial  electron  diffusion  across  a superimposed 
magnetic  field.  This  is  the  first  such  determination  in  a 
large-facility  BPD  and  points  to  a process  substantially 
faster  than  cross-field  collisional  diffusion.  The  reason 
for  this  is  identified  with  the  presence  of  large  amplitude 
plasma  waves  and  their  associated  electric  fields. 

While  enhanced  cross-field  diffusion  appears  to  be  a 
characteristic  of  many  turbulent  plasma  environments,  the 
observations  reported  here  are  particularly  unique  in  that 

3 

the  applied  magnetic  field  is  at  best  10  times  weaker 
than  in  any  other  turbulent  diffusion  investigation.  This 
points  to  the  importance  of  an  unequivocal  determination 
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of  the  turbulence  spectra  and  the  associated  instability 
processes.  This  work  is  currently  underway  at  NRL. 
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